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The Euler and Navier—Stokes equations revisited 

Peter Stubb(0 


The present paper is motivated by recent mathematical work on the incompressible Euler and 
Navier-Stokes equations, partly having physically problematic results and unrealistic expecta¬ 
tions. The Euler and Navier-Stokes equations are rederived here from the roots, starting at 
the kinetic equation for the distribution function in phase space. The derivation shows that the 
Euler and Navier-Stokes equations are valid only if the fluid under consideration is an ideal gas, 
and if deviations from equilibrium are small in a defined sense, thereby excluding fully nonlinear 
solutions. Furthermore, the derivation shows that the Euler and Navier-Stokes equations are 
unseparably coupled with an appertaining equation for the temperature, whereby, in conjunction 
with the continuity equation, a closed system of transport equations is set up which leaves no 
room for any additional equation, with the consequence that the frequently used incompressibil¬ 
ity condition V • v = 0 can, at best, be applied to simplify these transport equations, but not to 
supersede any of them. 


1. Introduction 

The Euler and Navier-Stokes equations belong to the oldest equations in physics. They have been 
the subjects of innumerable publications and are widespread over the textbook literature. It may 
appear, therefore, that nothing new could be written about these equations. 

The reason for a renewed interest in these equations has come from the mathematical side, primarily 
with regard to the Euler equation. A central question in these mathematical studies has been 
whether the energy of the system is conserved by the soltutions. Scheffer [1] and Shnirelman [2] 
showed the existence of solutions which seem to indicate that kinetic energy is created out of nothing. 
Other investigations (e.g., Constantin et al. [3], Cheskidov et al. [4], De Lellis and Szekelyhidi [5]), 
following up Onsager’s conjecture, deal with just the opposite, the dissipation of energy without 
having a loss term in the equations. 

From a physical point of view these results are unacceptable. How should simple hydrodynamics 
have the potential to break a most fundamental physical principle, the conservation of energy in a 
closed system? But then, how could it be that mathematically correct solutions of physical equations 
should be physically invalid? The origin will have to be found in the equations, and in a possible 
violation of the applicability limits of these equations. 

It is tacitly assumed in these mathematical treatments that the equations used are applicable to 
arbitrarily large perturbations, and further that the fluid is incompressible, to the extent that the 
incompressibility condition V ■ v = 0 is used as if it had the rank of a transport equation, qualified 
to be used in place of the closing equation for the pressure. The equations used in the quoted papers 
are 


— + (v • V)v = -Vp (1) 

V • v = 0 (2) 

where p = p/Nm , p is the pressure, N the particle number density (assumed to be constant), and 
m the particle mass. 
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The Navier-Stokes equation, on the other hand, has received renewed interest through the Mil¬ 
lennium Prize endowed by the Clay Mathematics Institute (Fefferman [6]). According to the Clay 
Institute’s judgement, solutions of 


— + (v • V)v = -Vp + n Av (3) 

in conjunction with eq. (2) belong to the seven most important open problems in mathematics. 
The prescribed task is to find a proof for the existence and smoothness of solutions of [(3),(2)], and 
the expectation is that thereby progress be made to unlock the secrets hidden in the Navier-Stokes 
equation. Realizing that the system [(3),(2)] is built on drastic physical simplifications, it is not 
obvious that a successful accomplishment of the mathematical task would have physical relevance, 
and indeed we will see that it has not. In (3), [i is the kinematic viscosity. 

It will be the purpose of the present paper to examine the physical consequences of the simplifying 
assumptions entering into the equations [(1),(2)] and [(3),(2)], respectively. 

2. Basic equations 

The mother equation for any transport equation is the continuity equation in phase space, describing 
the distribution function /(u, r ,t), normalized by f f d 3 u = N. Considering a one-constituent 
neutral fluid without external source or loss terms, we have 

^ + (u-v)-V/+ F ■Vuf= 5 [ t ( 4 ) 

where 5f/5t symbolically denotes the temporal change of / due to the action of collisions, and F 
is a force acting on the particles individually. In the given context, F represents the intermolecular 
force. The kinetic equation (4) relates to point particles, and to ideal gases if F is set to zero. 

It will be assumed that collisions are sufficiently soft to be non-ionizing. Another assumption will 
be that the particles have no internal degrees of freedom. The inclusion of internal excitations 
would mean that separate equations would have to be formulated for every excitation level, with 
terms describing the transitions between them. Such a treatment would be beyond the scope of the 
present work, and so we have to take it as given that the particles have only their three translatorial 
degrees of freedom which explains the occurrence of the factors 3/2 or 2/3 in some of the subsequent 
relations. 

Eq. (4) will be used to derive a system of transport equations for the moments of the distribution 
function. Before proceeding, we will distinguish two groups of moments, those originating from 
the distribution function in phase space, /, and those originating from the distribution function in 
velocity space, / = //IV. The macroscopic velocity v and the kinetic temperature T belong to the 
second group, 


v = / /ud 3 u 


T= — / /(u-v) 2 d 3 u 
3 K J y ' 


(5) 

( 6 ) 


where K is Boltzmann’s constant. A complete set of moments of / can be defined by 
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with 


( Fj = J /( Uj r5 *) ( u ) d 3 u (7) 

^,(*j,fe)(u) = m (u x - V x ) l (u y - Vy) J (u z - ^) fc (8) 

where i, j, k are non-negative integers, and n = i+j + k is the order of the moment. Multiplication 
of the terms in (4) by w j^] 1 subsequent intergration over velocity space yields, after some 

elementary mathematical steps, 


+ 

+ 


f .dv x ,dv v , dv z \ 

/ j ^(*- 1 d+ 1 ,fc) j M (i+i,j-i,k) ^y_ j j^(i,j-i,k+i) ® v y _|_ 

\ <9y cte 9s: 

ay ox oz J 

(- 

( dx dy 


f dM (* +1 J,fe ) 

."i 4” r\ 4~ 



(9) 


whereby the entirety of transport equations is concentrated in one equation. For the conversion of 
the last left-hand side term in (4) into a moment term it has been necessasry to postulate that the 
decrease of / is stronger than the increase of as u goes to infinity. 

The terms in (9) have been ordered by setting curly brackets. The first pair of brackets contains the 
wanted moment of order n, the second moments of order n— 1, the third other moments of order n, 
and the fourth moments of order n + 1. Eq. (9) shows that any wanted moment M i s coupled 
with a finite number of moments of same or lower order, and an infinite number of moments of 
higher order. This reminds us of the well-known, but often forgotten fact that transport equations 
do not exist as isolated entities, but only as integral parts within a system of transport equations 
which have to be extracted from the infinite system on the basis of precisely stated truncation 
conditions. 

Before making use of (9), we identify the moments of orders n = 0 to 3. 

For n = 0: 


Me 0 ’ 0 ’ 0 ) = Nm 

For n = 1: 

M mo) = M (0,1,0) = M (0,0,1) = o 


( 10 ) 


( 11 ) 


For n = 2: 


m ( 2 ’°’ 0) = fc 

, jVf( 0,2 ’ 0 ) = pyy , 

M(°-o.2) = p „ 

(12) 

M {;i ’ 1 ’ 0) = P X y = Py X , 

M ( I’O’l) = Pxz = Pzx 

, M (0 ’ ld) = Py Z = P Z y 

(13) 

T — 1 

3 KN 

(Ml 2 ’ 0 ’ 0 ) + M (0,2 ’ 0) + 

M(°’ 0 ’ 2 )) 

(14) 
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For n = 3: 


q x = ^(M^+M^+M^) 

(15) 

q y = ^(M^+M^+M^) 

(16) 

q z = l{M^+M^+M^) 

(17) 


These are the three components of the heat flux vector q. 


3. Transport equations for the fluid variables N, v and T 

Eq. (9), together with the moment definitions (10) to (17), will now be used to extract equations 
for N, v and T, neglecting intermolecular forces. We obtain: 

For n = 0: 


dN 

dt 


—IV V • v 


(18) 


For n = 1: 


d i = -jk lv(NKT)+v -^° 1 (19) 

where (p)° is the traceless part of the pressure tensor p, with elements defined by (12) and (13). 
The term in the first round brackets in (19) corresponds to the kinetic pressure 


p = NKT 


( 20 ) 


of an ideal gas. 
For n = 2: 


dT 

dt 


- 2 TV ■ v- 
3 


2 

3KN 


[(p)°: Vv + V-q] 


( 21 ) 


where Vv denotes the dyadic product. 

The reason why we have no collision terms in the transport equations for N, v and T is because 
we have assumed that collisions are non-ionizing ( SN/dt = 0), and because there exists no second 
constituent with which momentum and energy could be exchanged (5v/5t = 0 and ST/dt = 0). 
However, the randomizing effect of collisions will play a dominant role in the subsequent treatment 
of (p)° and q. 

There are two reasons why the equations (18) to (21) are not applicable to fluids in general, but 
only to ideal gases. The first is that intermolecular forces have been ignored. The second is that 
the whole formulation is built on the idealization that the particles are infinitesimally small, so that 
a given volume is able to accomodate an infinite number of particles. This becomes apparent from 
the continuity equation (18) which allows N to become arbitrarily large. 

In order to get an indication of the applicability limits of eqs. (18) to (21), we make recourse to van 
der Waal’s equation of state, written here in the form 
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Pt = Pk - Pi with p k 


NKT 
1-Vo N 


( 22 ) 


and pi = —N z 


where pt is the total pressure, p k the kinetic prassure, and pi the internal (or cohesive) pressure. 
The coefficient a, for various gases and liquids, is a tabulated quantity, A is the Avogadro number, 
and Vo is a measure of the volume of a single particle. 

We see that the numerator in the expression for p k represents the kinetic pressure of an ideal gas, 
whereas the denominator accounts for the fact that the particles have a finite volume, represented 
by Vo- Obviously, for having an ideal gas, the denominator must be close to 1, and thus NVo <C 1. 
The corresponding smallness condition for the internal pressure is ( a/A 2 )N 2 <C NKT. So we find 
that the equations (18) to (22) will be applicable if the two conditions 


IV < 


1 

% 


and 




KTA 2 

a 


(23) 


are satisfied. 

The above transport equations (19) and (21) are open-ended since the moments (p)° and q establish 
a link to an infinite number of higher order moments. It will be necessary, therefore, to find criteria 
allowing to express (p)° and q in terms of N. v and T , in order to obtain a closed system of transport 
equations with no other variables than these. This task will be addressed in sections 4 and 5 below. 


4. Approximation 1: Fully isotropized system 

If the dominant role of randomizing collisions goes to the extreme that the system is completely 
isotropized, the (p)° and q terms disappear, and we have 


dv 

dt 


1 

Nm 


Vp = -V 2 


VT VN 

~Y + ~W 


(24) 


dT 

dt 


- 2 rv-v 

3 


(25) 


where V is the thermal velocity ( V 2 = KT/m ). The first part of eq. (24) is the Euler equation, 
however with the restriction that p here is not the pressure in a general sense, but the kinetic 
pressure given by (20). The equations (18), (24), (25) form a closed set, leaving no room for any 
other equation for the fluid variables N, v and T. 

The derivation shows that the Euler equation has a very narrow range of applicability: Firstly, the 
fluid must be treatable as an ideal gas, in the sense of conditions (23), and secondly, deviations of 
the pressure tensor p from NKT U (U = unit tensor) must be exceedingly small (c.f. conditions (26) 
below). 


5. Approximation 2: Nearly isotropized system 

With reference to pressure tensor (p)°, complete isotropization means that p VJ = 0 (i ^ j) and 
pa = p. The plausible next step then is to go from completely isotropized to almost isotropized by 
demanding 


\pij | < Pn and | pa - p\ < p 


(26) 
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where p here stands as abbreviation for p = (1/3 ){p xx + Vyy + Pzz)- The same conditions apply to 
the Euler equation, but in a more rigorous way, which might be indicated by using the symbol 
instead of <C. 

In order to obtain local transport coefficients, i.e., transport coefficients which are determined 
entirely by the properties of the substance, independent of the temporal and spatial development 
of the process, it has to be demanded that the time between collisions, r c , is the shortest time, and 
the travel distance between collisions, l c , the shortest length in the system, 


t c <C t and l c <C l (27) 

where r is a characteristic time and l a characteristic length of the process. These conditions exclude 
extremely rapid and/or extremely short-scaled processes. 

We consider the pressure tensor elements p xx (i = 2,j = k = 0) and p xy (i = j = 1, k = 0). Since 
Sp/St = 0, the collision term 5p xx /5t in (9) can be altered into 5(p xx — p)/5t , and from the first line 
in (9), making use of (27), we obtain 

5{p xx -p) = d(NKT) + NKT 
St dt 

Next, expressing dN/dt by the continuity equation (18), and dT/dt. by the leading term in (21) 
( dT/dt = —(2/3)TV • v), we obtain 

=nkt U3v ! l _2 \ 

dt V dx 3 J v ’ 

From the first and third brackets in (9): 

= + <30) 

Corresponding expressions are obtained for the other elements of the pressure tensor, and the results 
can be written in compact form as 


V-v + 2 


dv x 

dx 


(28) 


*(P)° 

St 


NKT 


Vv + (Vv)‘ 


l (V • v)U 


(31) 


where (Vv)* denotes the transposed dyadic. 

The derivation of a corresponding expression for Sq^/St is a rather intricate matter since it will be 
necessary to go up to moments of order n = 4. For these we assume 


WoM k) \ « and - M&W I « 


(32) 


where M^en* means that all indices i, j, k are even, that not all indices are even, and Mmaxw 

that Meltin' 1 is approximated by using a local Maxwell distribution, displaced by the macroscopic 
velocity v. The background for these assumptions is that odd moments occur only in the form of a 


deviation from equilibrium, whereas even moments also exist in equilibrium. 
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We consider the x—component of q. From (9) with (15), (26), (27), (32) and F x = 0, adding the 
three equations for M^ 3,0,0 \ M^ 1,2,0 ^ and Af( 1,0,2 ), we obtain 

$Qx = 5 (fox 1 d_ ( m (4,0,0) , M (2,2,0) , M (2,0,2)\ (ooN 

5t 2 y dt 2 Qx \ maxw maxw 1 - LV - L maxw J 

With the useful relation 

= (* - 1)!! (j - 1)!! (k - 1)!! NKT V n ~ 2 (34) 

(where jbb.k) = 1 if z, j, k are all even ; S^ ,k ^ = 0 otherwise), the terms in brackets in (33) add up 
to 5NK 2 T 2 /m. Then, approximating dv x /dt by the leading term in (19), — (1/Nm) d(NKT)/dx, 
and generalizing from 5q x /5t to 5q/5t, we obtain 

= ^NKV 2 VT (35) 

ot 2 

To finalize the derivation, we need a relation between S(p)°/5t and (p)°, and between 5<\/8t and q. 
For this, we use the 10-moment relaxation model of Stubbe [7] and have 

ML = -2^ r \p)° (36) 

where the collision frequency is a measure of the strength of the randomzing action of collisions, 
quantitatively defined by eqs. (28c) and (15) of [7]. From (31) and (36) 


where 


(p)° = -v 


Vv + (Vv)‘ 


(V • v) U 


NKT 

71 ~ 2 z» 


is the dynamic viscosity. In a corresponding way, q is obtained as 


(37) 

(38) 


where 


q = — kV T 

_ 5 NK V' 2 
" “ 4 1A r ) 


(39) 

(40) 


is the heat conductivity. Since is proportional to N, both ?? and k are independent of N. The 
viscosity given by (38) agrees fully with the corresponding result provided by the first approximation 
of the Chapman-Enskog theory (see Chapman and Cowling [8]), whereas in eq. (40) the factor 5/4 
should be replaced by 15/8. 

We have now reached our goal to obtain a closed system of equations for the fluid variables N, 
v and T: The number density N is determined by (18), the macroscopic velocity v by (19) in 
conjunction with (37), and the kinetic temperature T by (21) in conjunction with (37) and (39). 
The applicability of this system of equations is bound to the conditions (27), (26) and (32). Whereas 
(27) is easily fulfilled under regular circumstances, (26) and (32) are strongly limiting since they 
have the character of a hidden linearization. 
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In order obtain the Navier-Stokes equation from (19) with (20) and (37), we have to assume that 
r] is constant (i.e., Vr/ sufficiently small). Hereby, (19) becomes 


dv 9 

' VT 

VN' 

71 

o 

— = -V 2 
dt 

~Y + 

N 

+ V* 

Nm 

Vv + (Vv) 4 — - (V • v) U 

o 


It is straightforward to convert (41) into the equivalent form 


(41) 


A = _y 2 

dt 


VT VN 

~Y + ~W 


+ 


Nm 


Av + - V(V • v) 

O 


(42) 


which is the Navier-Stokes equation, except that usually the first right-hand side term is written 
as — (1 /Nm)Vp, enabled by (20). 


Similarly, treating k as a constant, (21) in conjunction with (37) and (39) is converted into 


dT 

dt 


2 

3 


TV-v + 


2 TJ 

3 ~KN 


Vv : [Vv + (Vv)*] 


^(V v) 2 


2 

+ x 


K 


3 KN 


AT 


(43) 


Just as in the case of the Euler system [(18),(19),(21)], the equations (41) (or (42)) and (43), in 
conjunction with the continuity equation (18), represent an unseparable set, tied together by joint 
closure prescriptions. The number of equations matches the number of unknowns, leaving no room 
for any additional equation for the fluid variables N , v and T. And just as in the case of the Euler 
system, the Navier-Stokes system [(18),(41),(43)] is valid for ideal gases only. 

The validity of the momentum equation (41) and the thermal equation (43) is strongly restricted 
by the conditions (26) and (32). The consequence of these conditions is that the viscosity and heat 
conduction terms have to be small in comparison with the leading terms, i.e., the first terms each 
on the right-hand sides of (41) and (43). This offers an easy way to check mathematical solutions 
of the system [(18),(41),(43)] for physical validity. 


6. Energy equations 

As a useful addition, we convert (41) and (43) into energy equations. The kinetic energy density of 
an ensemble of particles due to its bulk motion is given by 


1 

wk = ~ m 


f v 2 d 3 u = -Nmv 2 
■' 2 


and the internal energy density by 


wi = - m I / (u — v) 2 d 3 u = | NKT = - p 
2 J 2 2 

From (19), multiplied by v, and (21), the energy equations 

+ v • (w K v) = - v [Vp +V- (p)° ] 


(44) 


(45) 


(46) 


dwi 


+ V • (wi v) 


[p V • v + (p)° : Vv] — V • q 
8 


dt 


(47) 



are obtained. With the identity v • [V • (p)°] + (p)° : Vv = V • [(p)° • v] , the sum of wk and wj 
w = wk + wi, is given by 


dw 

— + V • (rev) = -V • (pv) -V - [(p)° • v] - V q 

implying that the total energy f w d 3 r inside a closed solid surface is constant. 
Insertion of (37) and (39) in (46) to (48) yields 


(48) 


dwx 

~df 


+ V • [wk v) = — V • Vp + 7] V • V 


Vv + (Vv)* - - (V • v) U 


(49) 


dwi 

dt 


+ V • (iui v) = -pV-v|i] 


Vv : [Vv + (Vv)*] - -(V-v ) 2 
3 


+ k AT 


(50) 


+ V • (w v) = —V • (p v) + 77 V • ( v 


Vv + (Vv)* - - (V ■ v) U 


+ k AT 


(51) 


Any mathematical solution of the system of transport equations has to satisfy the energy equation 
(51) in order to be regarded as physically valid. 


7. Incompressible flow within the frame of the Navier—Stokes equation 

As we have seen, the Navier-Stokes equation does not exist as an isolated equation, but only as 
a part of a system of equations, based on well-defined joint applicability conditions, one of these 
being that the fluid under consideration is an ideal gas. Since an ideal gas is the least incompressible 
among all conceivable fluids, it appears pointless do devote much attention to this special case which, 
in the end, may have no relation to any physical reality. 

The incompressibility condition V • v = 0 is commonly used in combination with ViV = 0, 


V • v = 0 and VN = 0 (52) 

implying that N(r, t ) is constant and the continuity equation thus not needed. 

The conditions (52) are frequently used, but seldom backed by a solid foundation. A typical example 
of a quasi-justification, taken from Landau and Lifshitz [9], reads as follows: “In a great many of 
cases of the flow of liquids, and also of gases, their density may be supposed invariable, i.e. constant 
throughout the volume of the fluid and throughout its motion. In other words, there is no noticeable 
compression or expansion of the fluid in such cases. We then speak of incompressible flow.” But 
why should one not as well write “In a great many of cases, the temperature may be supposed 
invariable. We then speak of isothermal flow”? Both statements appear equally arbitrary. 

Despite these reservations, we will examine to what extent (52) will be able to simplify matters. In 
the literature, the conditions (52) are commonly used twofold, firstly to simplify the Navier-Stokes 
equation, and secondly to discard the temperature equation (43) and use V • v = 0 instead. The 
system [(3),(2)] is an example: The Navier-Stokes equation there is written in the simplified form 
(3), and eq. (2) stands where an equation for p should stand. 

The common way to employ the condition V • v = 0 in the Navier-Stokes equation (42) is to ignore 
the V(V • v) term, 
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( 53 ) 


dv 2 VT rj 

— = - V -1-— Av 

dt T Nm 

Realizing that Av = V(V • v) — V x (V x v), a more appropriate formulation would be 


— = -V 2 — 
dt T 


^Vx(Vxv) 

Nm 


A third approximation is obtained from (41), 


dv 

dt 



V 

Nm . 


V- [Vv + fVvf] 


(54) 


(55) 


In order to decide which of these three simplified momentum equations should be used, we go back 
to the origin, eq. (9). We see that the only occurrence of V • v there is in the first line. Setting this 
term to zero and following the derivation through to the Navier-Stokes approximation, we arrive 
at (55) which, therefore, is to be regarded as the appropriate one among the three versions. 

The thermal equation (43) with V • v = 0 becomes 


dT 

dt 


V 


3 KN 


Vv : [Vv + (Vv) 


2 

+ x 


K 


3 KN 


AT 


(56) 


The simplified equations (55) and (56) form a closed set, with N treated as an externally given 
constant. This is as far as one can get by using the assumptions (52) for a simplification of the 
momentum and thermal equation. 

The next step in the chain of approximations will be to discard (56) and use (2) instead. To see 
what the physical implications of this procedure will be, we go back to the thermal energy equation 
(47), write out the terms in squared brackets in component form, and define: 


WlzzP = -pV • v - (p xx 


dv x 


dv y dv z 


(57) 


( dv x 

WK^I = ~Pxy 


dVy\ 

dx ) 


Pyz 


f dvy dv z \ 

\ dz dy ) 


Pxz 


f dvz dv x \ 
\ dx dz ) 


(58) 


wi^i = k AT 


(59) 


Thereby, eq. (47) is brought into the concise form 


8wj 

~dt 


+ V • (wi v) 


Wl-^P + WK^I + Wl^I 


(60) 


The meaning of the energy exchange terms wi^p and is obvious: wp^p describes the mutual 

conversion of potential and internal energy, due to compression or dilatation, and wj^j describes 
the local redistribution of internal energy due to conduction of heat. With this given, the remaining 
term in (60), iuk^i, is left to describe the one-sided conversion of kinetic into internal energy due 
to internal friction. This interpretation is supported by the fact that vjk^i turns out to be an 
unconditionally positive quantity, as seen by insertion of (37) in (58): 
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W K ^I=V 


(61) 


(dv x dv v ^\ f dv v dv z ^\ (dv z dv x ^\ 

(l^ + fc) + (& + ^) + (ar + 97) 

The physical function of this term is to irreversibly transfer energy from the kinetic into the thermal 
energy reservoir. In a complete physical system, the gain term wk^i in (60) will be balanced by 
a corresponding loss term contained in the kinetic energy equation, and both will add up to zero. 
However, the system [(55),(2)] is not complete: The thermal energy equation is removed and replaced 
by the non-thermal equation V-v = 0, with the consequence that kinetic energy will be lost without 
appearing elsewhere. The total energy of the system will not be conserved. 

While wk->i is not affected by the conditions (52), wp^.p is. Insertion of (37) in (57) with (52) 
yields 


wi^p 



(62) 


The conditions (52) are meant to exclude compressions or dilatations, implying that wp-^p should 
vanish, yet we see from (62) that wj^p in fact does not fully disappear. The remaining part (62), 
being entirely positive, acts as if compressions are possible, but not dilatations. This is another 
unacceptable consequence. 

A third inconsistency is that use of the system [(55),(2)] includes the possibility of a local accumu¬ 
lation or depletion of thermal energy, but contains no mechanism for a redistribution by means of 
heat conduction. 

We conclude that the replacement of the thermal equation (56) by the incompressibility condition 
V • v = 0 has unacceptable consequences, even in the case that this condition as such might be 
justifiable. All that this condition can possibly do is to bring the full system of transport equations 
into the simplified form [(55),(56)]. As a consequence, equations [(3),(2)] will not be able to describe 
a real physical situation. Their solutions will be of sole mathematical interest. 


8. Incompressible flow within the frame of the Euler equation 

Repeating from section 4, the Euler system is given by 


<9v 

dt 


+ ( v • V)v = -V 2 


VT VN 
~T ~ + ~N~ 


1 

Nm 


Vp 


dT 

~dt 


+ v-VT 


--TV-v 
3 


with p = NKT and N to be determined by the continuity equation (18). 
Simplifying [(24),(25)] by means of the incompressibility conditions (52), we have 


dv 

dt 


+ (v • V)v = -V 2 


VT 

~T~ 


dT 

~dt 


+ v • VT = 0 


(24) 

(25) 


(63) 

(64) 
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These equations are now independent of N, and the temperature gradient is the only driving agent. 
Again, we have a closed set of equations, describing the fluid variables v and T, and again there is 
no place for any additional equation. 

The fact that these equations are based on the assumption V • v = 0 does not imply that their 
solutions will actually confirm the validity of this assumption. Therefore, if {vi(r, t), T\(r, £)} is a 
mathematical solution of [(63),(64)], it will be necessary to show a posteriori that V • vi ~ 0 for all 
r and t.. 

In this situation it may be tempting to enforce incompressibility by using eq. (2), V • v = 0, in 
place of (64). Although this is common practice, as exemplified by [(1),(2)], it does not make 
things easier, because now it has to be shown that any solution {v 2 (r, t), ^(r, £)} of [(63),(2)] obeys 
the temperature equation (64). This equation does not lose its existence only because the fluid is 
assumed to be incompressible. 

In the end we see that the full Euler system [(24),(25)] is easier to handle than its simplified versions 
[(63),(64)] or [(63),(2)]. Incompressibility is not a property that can be enforced; if existing, it will 
have to result self-consistently from the full equations. 

The arguments presented here apply correspondingly to the Navier-Stokes system, in addition to 
those given in section 7. 


9. Summary 

It has been the purpose of the present paper to work out the applicability limits of the Euler and 
Navier-Stokes equations, and to examine the consequences of the frequently used incompressibility 
assumption. To this end, the Euler and Navier-Stokes equations have been rederived, starting at 
the kinetic equation for the distribution function in phase space. 

It follows from this derivation that the applicability of the Euler and Navier-Stokes equations is 
bound to the following conditions: 

(a) The localization condition (27), securing that the transport coefficients are matter constants, 
and excluding extremely rapid and/or short-scaled processes. 

(b) The isotropization conditions (26) and (32). These have the character of a hidden linearization, 
implying that only solutions desribing small deviations from equilibrium are physically valid. 

(c) The rarefaction conditions (23). These exclude the applicability of the Euler and Navier-Stokes 
equations to dense gases and liquids. For extremely rarefied gases, a conflict may arise between 
conditions (a) and (c), with the consequence that fluid theory becomes inapplicable altogether. 

(d) The energy condition, consisting in the requirement that the energy equation (51) be satisfied, 
corresponding to the conservation of the total energy in a closed system. If mathematical solu¬ 
tions are obtained that do not meet this requirement, it will be pointless to search for physical 
explanations, there are none. 

The conditions (a) to (d) set a dividing line between physically valid and invalid solutions, and this 
dividing line must be observed when hydrodynamic equations are used as material for mathematical 
studies. 

Further, it follows that the system of equations comprising the Euler or Navier-Stokes equation 
forms a complete set for the fluid variables N, v and T, leaving no place for any other equation, 
like the relation V • v = 0 used to enforce incompressibility. Consequently, the equations [(1),(2)] 
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or [(3),(2)] are without physical foundation and must be left to purely mathematical interest. 
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